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Abstract. Single crystals of the heavy- fermion compound YbAgGe have been 
studied by neutron scattering. The magnetic ordering occurring below T% ~ 0.5 
K is characterized by a commensurate propagation vector k= (1/3,0, 1/3) and the 
moments in the basal plane of the hexagonal structure. The dynamic magnetic 
susceptibility is dominated by quasielastic spin fluctuations with a characteristic energy 
r q (T) of the order of 1 meV. The spins fluctuate predominantly in the basal plane. 
No spin-wave excitations are observed in the magnetically ordered phase. Below 
the Kondo temperature, Tk ~ 20 K, T q (T) shows a strong q dependence for wave 
vectors along the c* direction, but is (/-independent in the basal plane. r q (T) shows 
initially a rapid increase with temperature T at the antiferromagnetic zone center, 
but follows a standard VT law for other q values and for T > Tk in general. 
These observations classify YbAgGe as a well-behaved heavy-fermion compound with 
a particular q-dependence of the antiferromagnetic spin fluctuations, possibly related 
to the geometrical frustration of the Yb 3+ ions. 
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1. Introduction 

YbAgGe is a recently discovered Yb-based moderate heavy-fermion (HF) compound 
with a Kondo temperature of T# ~ 20 K and a linear coefficient of the specific heat, 
7 = Cp/T, of 570 mJK~ 2 mole -1 at T = 1.6 K and 150 mJK _2 mole _1 in the zero- 
temperature limit [HEJIH]- At low temperatures, YbAgGe shows two magnetic phase 
transitions, a first-order transition at T\ ~ 0.65 K and a second-order transition at 
T 2 = 0.8-1.0 K [H Hj. The small entropy at low temperatures, only 5% of i?m2 at 
T = 1 K, suggests a small-moment ordering pQ. Both magnetic phases disappear with 
the application of a magnetic field of the order of 1-3 and 4-8 T, depending on the 
field orientation [21 H|- At the same time, non- Fermi liquid behavior is observed, with a 
logarithmic term in the specific heat, C p /T oc — InT, and a linear resistivity, p — Po ^ T, 
for fields of the order of 8-12 T E] • YbAgGe provides thus an excellent opportunity 
to study a field-induced non-Fermi liquid of a stoichiometric compound in the vicinity 
of a quantum phase transition (QPT). We have performed neutron scattering studies of 
single-crystalline samples in the zero-field heavy-fermion regime, which is characterized 
by both the Kondo effect and geometrical frustration. The measured dynamic magnetic 
susceptibility has an energy and temperature dependence that is typical for a heavy- 
fermion compound, but the wave vector dependence is unusual. 

YbAgGe crystallizes in the hexagonal ZrNiAl-type crystal structure with the 
non-centrosymmetric space group P62m and room-temperature lattice parameters of 
a = 7.05 and c = 4.14 A (see figure jOj. The Yb 3+ ions have 4/ 13 electronic 
configuration, which is the hole counterpart to trivalent Ce (4/ 1 ) ions. The three Yb 
ions per unit cell occupy the same equivalent site in the lattice and lie on a Kagome- 
like triangular lattice in the basal plane, which for antiferromagnetic coupling leads 
to geometrical frustration. The eight-fold J = 7/2 multiplet of the Yb 3+ ion is split 
into 4 doublets by the crystal electric field. Specific-heat data show a peak at T = 60 
K, which is attributed to the Schottky contribution from the first excited doublet at 
an energy of E = 9.4 meV jHj- However, the peak height is two times higher than 
expected for a doublet-doublet excitation, implying that the third and fourth doublets 
also need to be taken into account. The magnetic susceptibility, x-, follows a Curie- 
Weiss law at high temperatures with an effective moment of p e ff = 4.4 \xb PHE], i-e. 
close to the value of the free Yb 3+ ion, p e ff = 4.5 and a low-temperature easy- 
plane anisotropy of Xab /Xc ~ 3 [TJ [21 E] . A small maximum at T = 4 K is seen in Xa 
PU El] , reflecting the onset of antiferromagnetic correlations. The corresponding crystal- 
field scheme, which has the a-axis as principal axis, has been confirmed by inelastic 
neutron scattering measurements on a polycrystalline sample, which shows a clear peak 
at E — 12 meV [7j. These measurements also revealed Kondo-type spin fluctuations 
with a characteristic energy of T = 0.9 meV at low T and a \/T behavior at higher T. 

We have performed neutron scattering measurements on large single crystals of 
YbAgGe. The sample growth and experimental conditions are discussed in section [2] 
SectionElpresents neutron scattering measurements of the magnetic structure of the low- 
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Figure 1. Crystallographic structure of YbAgGe. 

T phase below T\. Inelastic neutron scattering measurements of the dynamic magnetic 
susceptibility are presented in section |U and the results are discussed and compared with 
other Af heavy-fermion materials in section 03 

2. Experimental 

Two single crystals, each of approximate weight of 2 g, were grown using a Ag-Ge rich 
self- flux method fl]. Smaller crystals using the same growth technique have been fully 
characterized by macroscopic measurements QjJ E] . Neutron scattering measurements 
were performed on the single-crystal white-beam time-of-flight diffractometer SXD at 
the ISIS spallation source and on the cold triple-axis spectrometer IN14 of the ILL high- 
flux reactor. In this paper, we label the coordinates in reciprocal space using hexagonal 
reduced lattice units (r.l.u.) with indices (hkl), where a point in reciprocal space is 
given by Q= ha* + A;b* + Ic* with a* = An/ (y^a) x and c* = (27r/c) z. 

On SXD, one single crystal was aligned with the c axis vertical and glued by black 
Stycast to a Cu support. The assembly was mounted in a 3 He refrigerator which was 
inserted in a standard helium cryostat. Both the cryostat and the 3 He insert were 
equipped with Al tails. Measurements were performed for three different orientations 
of the crystal (rotated around the c axis) at two temperatures: T = 0.32 K and T = 2 
K. All integer reflections could be indexed in the P62m space group. The magnetic 
scattering was obtained by subtracting the T = 2 K data from the low-T data, which 
removed nuclear Bragg reflections from the sample and powder lines from the cryostat. 

On IN14, two single crystals of YbAgGe were aligned together and mounted on a Cu 
support using Cu straps and Cu mesh for thermal contact. The assembly was mounted 
in a dilution fridge with (hOl) as horizontal scattering plane. IN14 was used in long- 
chair configuration with a vertically focusing PG(002) monochromator, a horizontally 
focusing PG(002) analyzer, and horizontal collimations of guide-60 '-open-open. Most 
measurements were performed at fixed final wave vectors of kf = 1.97 and kf = 1.30 
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Figure 2. Color map of the low-T magnetic scattering in YbAgGe, taken as the 
difference between data at T = 0.32 and 2.0 K, in the (hk) scattering plane for I = f /3. 
The six-fold symmetry of the magnetic satellites characterized by a propagation vector 
of k=(l/3, 0,1/3) is easily seen. 



A -1 , with energy resolutions for elastic incoherent scattering of 0.442 (2) and 0.117 (1) 
meV, respectively. Higher-order harmonics were removed from the scattered beam by 
pyrolytic graphite or liquid-nitrogen cooled Be filters. The intensity of the magnetic 
scattering so obtained is proportional to the dynamic structure factor S(q,u), 

S(q,w) = z \ Tps X"(q>^)> (!) 

1 — exp(— nu/k B T) 

where the dynamic magnetic susceptibility, x"(q, u;), is obtained from 5(q, to) after 
subtracting a constant background term, determined at negative energy transfers at 
low temperatures or at positive energy transfers with the analyzer turned off from the 
Bragg condition by 10 degrees. Measurement were performed at temperatures (of the 
mixing chamber) between T = 0.06 and 50 K in zero magnetic field. At the lowest 
temperatures, the sample could have been slightly warmer (~ 0.1 K) than the mixing 
chamber due to beam heating. 



3. Magnetic order 

The magnetic scattering in the low-temperature phase at T = 0.32 K (after subtraction 
of the T = 2 K data) measured on SXD is shown in figure 121 Non-integer superstructure 
peaks due to magnetic order are clearly observed. They can be indexed with a 
propagation vector of k=(l/3, 0,1/3). The intensity of the superstructure peaks 
decreases with increasing Q, as expected from the magnetic form factor. The intensity 
is weak for small I indices, which suggests that the moments are predominantly in the 
basal plane of the hexagonal structure, in agreement with the easy-plane anisotropy 
observed in the magnetic susceptibility measurements. 

Scans of the magnetic Bragg peaks along the h and / directions were performed 
below and above the magnetic ordering temperature T\ on IN14 at kf = 1.3 A -1 
(see figure EJ). These scans confirm the propagation vector of k=(l/3, 0,1/3) observed 
on SXD. The magnetic order is long-range, as seen from the resolution-limited width 
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Figure 3. Scans at zero energy transfer along the (a) h and (b) I direction of the 
Q=(l/3, 0,2/3) magnetic Bragg peak of YbAgGe measured on IN14 using kt = 1.3 
A -1 . Closed and open circles correspond to temperatures of 0.1 and 5.0 K, respectively. 
The lines are Gaussian fits. 
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Figure 4. Temperature dependence of the magnetic Bragg peak intensity at Q= 
(1/3,0,2/3) of YbAgGe measured on IN14 using k f = 1.3 A" 1 . The line is a guide to 
the eye. 

of the peaks. The temperature dependence of the magnetic Bragg peak intensity is 
shown in figure EJ The transition temperature determined from these measurements is 
T\ = 0.48 (2) K, which is similar to that from bulk measurements. The rapid drop of 
the magnetic intensity near T\ is reminiscent of a first-order phase transition. Above 
Ti, no magnetic super-lattice peaks were observed. Further measurements are planned 
to investigate if the phase between T\ ~ 0.5 K and T 2 ~ 0.9 K corresponds to a 
magnetically ordered phase. A full determination of the magnetic structures below T\ 
is also in progress. 

4. Dynamic magnetic susceptibility 

Our inelastic neutron scattering (INS) measurements on single crystalline YbAgGe show 
the existence of strong magnetic quasielastic scattering arising from spin fluctuations. 
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Figure 5. Energy scans of the quasielastic magnetic scattering at two equivalent 
antiferromagnetic positions, measured on IN14 at T < 0.1 K using kf = 1.97 A -1 . 
The lines are fits to a quasielastic Lorentzian, Eqs. lllll'l) . plus a Gaussian describing 
the elastic scattering. 



The dynamic magnetic susceptibility x"(q, u;) is quite well described by a quasielastic 
Lorentzian, 



with a characteristic energy of T q ~ 1 meV. This is in agreement with neutron scattering 
measurements on a polycrystalline sample by Matsumura et al [7j. Matsumura et 
al also observed a crystal-field excitation at approximately 12 meV, which is outside 
the energy window of the present measurements. The intensity of the spin-fluctuation 
scattering at the antiferromagnetic position Q=(l/3, 0,4/3) is nearly twice as big as 
at the equivalent position Q=(4/3, 0,1/3) (Figure EJ). Since INS probes magnetic 
fluctuations perpendicular to the total wave vector transfer Q, this suggests that the 
spins fluctuate predominantly in the basal plane of the hexagonal structure. This is in 
agreement with the easy-plane anisotropy in bulk magnetic susceptibility measurements. 

As we will show in detail below (section I4.2J1 . the quasielastic magnetic scattering is 
independent of temperature between T = 0.1 and 0.5 K, i.e. the quasielastic scattering 
exists in the low-temperature antiferromagnetically ordered phase. This may be a 
feature common to 4/ heavy-fermion materials close to a quantum phase transition, 
as e.g. quasielastic magnetic scattering has also been observed in CePd 2 Si 2 |E| and 
Celn 3 jOj in the ordered state. In these compounds, the spin fluctuations coexist with 
spin waves. In YbAgGe, we have not observed any spin waves. This may be because of 
the weak magnitude of the ordered moment. 

4-1- Q dependence 

Figure |H1 shows Q scans in two different directions taken at an energy of E = 1 meV, 
where the intensity of the spin-fluctuation scattering is maximal. The intensity at 
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Figure 6. Scans in wave vector Q along the (a) h and (b) I direction of the spin- 
fluctuation scattering at an energy transfer of E = 1 meV, measured on IN14 using 
kf = 1.30 A -1 . Closed and open circles correspond to T < 0.1 and T = 27 K, 
respectively. The lines are guides to the eye. The inset shows the scan directions in 
the (hOl) plane. 




Figure 7. Energy scans of the quasielastic magnetic scattering at low temperatures, 
T < 0.1 K, for different q values along the (a) h and (b) I direction, measured on IN14 
using kf — 1.30 A" 1 . The lines are fits to a quasielastic Lorentzian, Eqs. (11 1211 . plus a 
Gaussian describing the elastic scattering. 
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Figure 8. Wave vector dependence of the static susceptibility x'q and the characteristic 
energy T q of the quasielastic magnetic scattering along the h (left part) and / (right 
part) direction at T < 0.1 K (closed circles) and T — 27 K (open circles). The lines 
are guides to the eye. 



low temperatures is clearly modulated in the I direction. This shows that the spin 
fluctuations are correlated in Q and hence does not correspond to localized (or single- 
site) fluctuations. Such correlated spin fluctuations have been observed in other heavy 
fermion systems. The modulation of the intensity is not simply a modulation in 
but corresponds also to a q dependence of the characteristic energy T q . This is clearly 
seen in energy scans at different Q values along I taken at low temperatures (figure 
Eb)- On the other hand, the Q scan taken in the h direction (cf. figure EJ) is nearly 
flat, and energy scans along the h direction (cf. figure [T^l) are nearly the same. The 
corresponding "dispersion" of the characteristic energy T q as well as the q-dependent 
susceptibility x' q , obtained from fits of Eqs. (HH21), are shown in figure |H1 Clearly, both 
quantities are strongly dependent on q along the / direction but not along the h direction. 
Such a behavior of T q have not been reported in other heavy-fermion compounds and 
will be discussed in sectional 

4-2. T dependence 

The temperature dependence of the dynamic magnetic susceptibility was studied for 
two wave vectors: at the antiferromagnetic zone center Q! = (l/3,0,4/3) = (0,0,l)4J<; and 
at the AFM zone boundary Q =(l/3,0,1). These data, shown in figure El were fitted 
by a quasielastic Lorentzian, Eq. ©. The extracted temperature dependence of the 
characteristic energy T q (T) and the static susceptibility x q (^) are shown in figure 
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Figure 9. Dynamic magnetic susceptibility for different temperatures at (a) the AFM 
zone center Q=(l/3, 0,4/3) and (b) the AFM zone boundary Q=(l/3,0,1), measured 
on IN 14 using kf = 1.30 A -1 . The lines are fits to a quasielastic Lorentzian, Eq. J3J). 

Hm The characteristic energy at the AFM zone boundary, Fq, is constant at low 
temperatures and follows the standard x/T relation for heavy-fermion materials at higher 
temperatures, i.e. 

r (r) = r (o) + q(t -t*)aVt (3) 

with r (0) = 1.87 (7) meV, A = 0.17 (2) meVK" 1 / 2 , and T* = 5.0 (2) K. Here, T* 
that occurs in the step function is close to the temperature where the bulk magnetic 
susceptibility Xa peaks j2]. The characteristic energy at the AFM zone center, Ti, 
increases much more rapidly at low temperatures and joins Fo at T « 15 K, i.e. in 
the vicinity of T# ~ 20 K. Above TV, it follows the same a/T relation as r . The 
q-dependent magnetic susceptibilities x' q also show a similar T dependence at the two 
wave vectors, with (very approximative) oc 1/r (see section |S] for a more detailed 
discussion). At temperatures above T K ~ 20 K, the magnetic scattering at the AFM 
zone center is very similar to that at the zone boundary, showing that the correlations 
in q vanishes at high T. This is also clearly seen in the Q-scans at E = 1 meV, where 
the intensity modulation has disappeared at T = 27 K (see Figure EJ). 

5. Discussion 

As a consequence of the large moment in the doublet crystal-field ground state, the 
intensity of the low-energy quasielastic magnetic scattering from spin fluctuations in 
YbAgGe is relatively strong, which allows a detailed analysis of the q dependence of the 
dynamic magnetic susceptibility, w )- We find that both %q and T q depend not only 
on the magnitude q — |q| of the wave vector but also on its direction (see figure|HJ). There 
are only few examples of such behavior in 4/ heavy fermion materials. In cerium HF 
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Figure 10. Temperature dependence of the static susceptibility X q (T) and the 
characteristic energy r q (T) of the quasielastic magnetic scattering at the AFM zone 
center Q=(l/3, 0,4/3) (solid circles) and at the AFM zone boundary Q=(l/3,0,1) (open 
circles). The dashed line shows a fit of Eq. © to the zone boundary energy r . 



compounds, a q-dependent x"(q, has been observed in single crystalline samples of 
CeCu 6 HDlini, CeCu 5 . 9 Au .i [12], CeRu 2 Si 2 EJ HSJ , (Ce .925La .o75)Ru 2 Si 2 PIE], 
and CeNi 2 Ge 2 [TH]. However, little is known about the q dependence of the characteristic 
energy T q in these systems, because of the limited statistics in the neutron scattering 
measurements. An explicit q (but not q) dependence of T q has only been reported in 
CeCu6 [Ej and CeRu 2 Si 2 [13]. CeNi 2 Ge 2 is an exception: the dynamic susceptibility has 
two components, and for the high-energy (E — 4 meV) response, T q is independent of q 
while x'q shows a quasi two-dimensional behavior |18j . Some limited information about 
the g-dependence of u ) have also been extracted from powder data (examples 
are CeRu 2 (Sii_ :r Ge :c ) 2 an< i Ce(Ruo. 2 4Fe .76)2Ge 2 20 ). but most work assume that 
x"(q, a;) is either q independent or report quantities that are averaged over the Brillouin 
zone. 

In spin- fluctuation theories [21], the presence of an exchange interaction, J(q), 
between localized moments described by the local dynamic susceptibility Xl{w) = 
Xl/ (I — iu/Tl), can be described by an RPA-like expression x -1 ^ w ) = Xl}^) ~ ^(q), 
which leads to a quasielastic Lorentzian [Eq. (J2J)] with x q = Xx/[1 — J(<1)x'l\ an d 
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Figure 11. The product %T of the static susceptibility Xq(^) an d the characteristic 
energy T q (T) of the quasielastic magnetic scattering, (a) as a function of q along 
h (left part) and I (right part), (b) %T as a function of temperature for q at the 
antifcrromagnetic zone center and at the zone boundary. 



T q = Tl[1 — J(q)x' L ]. It then follows that the product XqF q = x'iXl becomes 
q independent, and also temperature independent (22J- Thus, in these models, the 
behavior of T q follows simply from that of x q - Experimentally, this relation has 
been approximately verified in CeCu 6 [TUj and CeRu 2 Si 2 [THJ ITS] . In YbAgGe at 
low temperatures, the product x'qFq is constant only in the h direction, where no 
q dependence of x"(q, to) is observed, while it varies with q in the "dispersive" I 
direction (see figure ITTk) . Also, the product x q (^)r q (T) is not constant as a function 
of temperature, cf. figure ITTb. 

As said, the most characteristic feature of the spin dynamics in YbAgGe is that 
T q depends on the direction of the wave vector q. The strong variation of T q and 
Xq for q along the hexagonal c axis is due to antiferromagnetic correlations. Similar 
antiferromagnetic correlations are also expected in the basal plane. However, neither T q 
nor x'q depend on q for q in the basal plane. A likely origin for this quasi-one-dimensional 
behavior is the geometrical frustration arising from the triangular coordination of the 
Yb ions, which lie on a quasi-Kagome lattice. Such frustration is expected to average 
out the correlations in the basal plane. However, it cannot be excluded that the quasi- 
one-dimensional behavior arises from Fermi surface effects. In addition to the particular 
q dependence of T q and x q , we also find that the spins fluctuate predominantly in the 
basal plane, due to the anisotropy of the crystal-field interactions. 

With increasing temperature, the characteristic energy of YbAgGe increases much 
faster at the magnetic zone center than at the magnetic zone boundary (see figure ITUj) . 
Above Tk ~ 20 K, the modulations in q disappear and all T q become identical. In 
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this temperature region, T q follows the yT dependence characteristic for heavy-fermion 
compounds. Similar behavior has been found in other heavy-fermion systems, such as 
CeCu 6 HI], CeRu 2 Si 2 [H], (Ce .925La .o75)Ru 2 Si 2 |T7] , and Ce(Ru . 2 4Feo.76) 2 Ge 2 [201, 
where T q= k at the critical wave vector k joins from below T q ^k at temperatures of the 
order of the Kondo temperature Tr. 

For a system close to a quantum critical point, the dynamic susceptibility is 
expected to follow so-called uj/T 13 scaling, 

x"(q,^)T a = f(u>/T% (4) 

where / is a scaling function, a ~ 1, and (3 = 1-1.5 [2H]- This implies that the only 
energy scale in the system is the temperature, and in particular, that T q (T) — ► as 
T — > 0. This is clearly not the case for YbAgGe in zero magnetic field, where T q is finite 
at all q in the limit T — ► 0. YbAgGe in zero field is thus a well-behaved heavy-fermion 
without u/T 13 scaling. 

The next step in our investigation of single- crystalline YbAgGe is to study by 
inelastic neutron scattering the evolution of the low-temperature dynamic magnetic 
susceptibility as the system is tuned by an applied magnetic field from the Fermi-liquid 
regime to the non- Fermi-liquid regime, where u/T 13 scaling might be observed. 
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